International Journal of Advances in Applied Sciences (IJAAS) 
Vol. 7, No. 1, March 2018, pp. 86~96 
ISSN: 2252-8814, DOI: 10.1159 1/ijaas.v7.11.pp86-96 o 86 


Experimental and Modeling Dynamic Study of the Indirect 
Solar Water Heater: Application to Rabat Morocco 


Ouhammou Badr’, Azeddine Frimane?, Aggour Mohammed’, Brahim Daouchi’, Abdellah Bah’, 
Halima Kazdaba’ 
1:234 Renewable Energy and Environment Laboratory, Department of Physics, Faculty of Sciences Ibn Tofail University, 
B.P 133, 14 000, Kenitra, Morocco 
ni Equipe de Recherche Thermique, Energie et Environnement (ERTEE), Laboratoire LM2PI ENSET, Universite 
Mohammed V, Rabat, Marocco 


Article Info ABSTRACT 

Article history: The Indirect Solar Water Heater System (SWHS) with Forced Circulation is 
, modeled by proposing a theoretical dynamic multi-node model. The SWHS, 

Received Jun 17, 2017 which works with a 1,91 m2 PFC and 300 L storage tank, and it is equipped 

Revised Jan 20, 2018 with available forced circulation scale system fitted with an automated sub- 

Accepted Feb 11, 2018 system that controlled hot water, is what the experimental setup consisted of. 


The system, which 100% heated water by only using solar energy. 
The experimental weather conditions are measured every one minute. 
Keywords: The experiments validation steps were performed for two periods, the first 
one concern the cloudy days in December, the second for the sunny days in 
May; the average deviations between the predicted and the experimental 
values is 2 %, 5 % for the water temperature output and for the useful energy 


Dynamic multi-node model 
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Solar collector application are 4 %, 9 % respectively for the both typical days, which is very satisfied. 

Thermal performance The thermal efficiency was determined experimentally and theoretically and 

Water heater shown to agree well with the EN12975 standard for the flow rate between 
0,02 kg/s and 0,2kg/s. 
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1. INTRODUCTION 

The most common use of thermal solar energy has been for water heating systems; this use has been 
commercialized in many countries in the world. Solar heating systems have many innovations and 
developments in the last years. Their applications have increased significantly, especially in countries with 
large solar potential such as our case Morocco. 

Thermal solar collector systems (water heating) can be used in two ways, by forced circulation 
systems or systems with natural circulation [1]. The modern research in thermal engineering is a combination 
of numerical and experimental aspects. In addition, many models and studies were carried out; for example, 
Duffie and Beckman [2] developed a model of design and prediction for the flat plate solar collector. In 
addition, Close et al [3], have developed a model named by In-node, also we have the model 2n-node and 
3n-node developed by [4] and [5] respectively, in the same direction the works of Isakson et al [6], De Ron et 
al[7], Huang et al [8], Schnieders et al [9], Fraisse et al [10], Cadaflach et al [11], and Molero Villar et al [12] 
have participated effectively in the development and modeling of this system. Furthermore, many norms 
describe the test procedures for the performances of thermal solar collectors have been generated. The best 
known are the [SO-9806-1 [13], the ASHRAE standard 93 [14] and EN12975 [15]. 

The simulation using the current model was shown to agree well with the experiment results. 
The objective of this work is to validate the theoretical model by performing a number of simple experiences 
for different weather conditions in purpose to evaluate the thermal performance of the SWH also, thinking of 
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using new alternatives for numerical simulation and use it after for another application as coupled it with 
others software (TRNSYS) in order to use the model in the building heating. 


2. MODELING OF THE SYSTEM 
2.1. The Heat Transfer Coefficients 

Before forming the model, we exposed the most used models in the estimation of the heat transfer 
coefficients. 

The transfer coefficient by convection between glass cover and ambient is calculated by the 
empirical relationship of Adams and Woertz [19]: 


h., =5,67 + 3,86V,, 


cv,out 


(1) 


The radiation coefficient between the glass cover and the sky is given by (Soteris Kalogirou 2014): 


h =¢,0(T, +T (T? +T) (2) 


r,g—sky 


Where the temperature of the sky is correlated by the following formula (Soteris Kalogirou 2014): 
Ty =0,0552 Ti; (3) 


The radiation coefficient between two infinite parallel surfaces, (1) and (j) [1]: 


LTE? -TP 4) 


y 
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The coefficient of conduction is given by the following equations: 
e; 
h cd = > K. (5 ) 


For the coefficient of convection between cover and absorber is calculated by the following empirical 
formula given by [1]. 


Nu K; 
hepe = ar (6) 
PTE D 
+ . 1,6 0,333 a 
vi Laadi = 1708 | hı _ 1708[sind.82)] |. Ra = + (7) 
Ra cos( p) Ra cos( p) 5830 
The Rayleigh value (Ra) is given by: 
Ra = =a, T (8) 


The forced convection coefficient absorber-fluid h, is calculated according to the empirical formulas 


Vif =P 
(equations 15, 16 and 17) given by the works of Haussen and Sider-Tate [19]. 
For Re < 2100: 
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Gz < 100 Nu = 3,66 + (9) 
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Gz > 100 Nu =1,86 Gz? + 0,87(1 + 0,015 Gz?) (10) 


For 2100< Re < 10000: 


2 1 


Nu =0,116(Re? 125) Pr? (1+——_) 
NL 


w| N 


(11) 


2.2. The Heat Transfer Coefficients 
As the other models, the present one is based on mass and energy conservation laws. The solar 

water heater is formed by many elements as shown in table.2 and in Figure. 1. 
Before starting, the following assumptions are made to simplify the analysis: 
The physical properties and thermal fluid data are based on the average temperature. 
Irradiation is uniform for all area of the collector. 
The ambient temperature is the same around the collector. 
The thermal exchange by convection between the insulation and ambient is negligible. 
Properties of materials are independent of temperature. 

Thus it is necessary to define step by step the singular heat flow equations in order to find the 
governing equations of the collector system. 


ono oR 


2.2.1. Heat Balance Of The Glass Cover 

The energy comes from the sun in the form of radiation and penetrates the aperture plane of the 
solar collector. However, a part of this radiation is reflected back to the sky, another component is absorbed 
by the glazing and the rest is transmitted through the glazing and reaches the absorber plate as short wave 
radiation. The equation 12 expressed the heat balance in this element. 


(m Cp), Ht “= sla, G) F h, „8-— ae Go 7 T,) + TE eas hy sai) os =l, )— (h, ,8— abs T hy ,g8— pT, Sd day 


2.2.2. Heat Balance of the Absorber 


The absorber receives energy coming from the sun through glass cover by radiation exchange with 
glass and at the same time transferred this energy into the fluid. Knowing that a part of energy transferred is 
lost because of the heat exchange with cover and with insulation by convection and conduction respectively 
as shown in the equation 13. 


T, = Pale u~ Di haa, p- Pee ins -T,)| a3) 
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dT, 
(mCp), —*=S|a,r,a,,G)-(h 
dt 
2.2.3. Heat Balance of the Fluid 


In this model, the exchanges existed only by forced convection between absorber and the water 
(working fluid) as it’s shown in the following equation: 


dT 
(OV Cp) p= ey, p-¢ ST =T) (14) 


2.2.4. Heat Balance of the Insulation 
The heat exchange concerning this element is those between the absorber and insulation by 


conduction (h ), also between insulation and ambient by radiation (h q ) as illustrate in the 


cd ,abs—ins r,ins— groun 


following equation: 


dT. 
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2.3. Performances of Solar Collector 

The thermal performance of solar collectors can be determined by the details analysis of the optical 
and thermal characteristics of the collector materials and collector design (theoretical modeling) or by 
experimental performance tested under different weather conditions (experimental validation). Solar 
collectors can be tested by two basic methods: under steady-state conditions or using a dynamic test 
procedure which adopted this work. 
The collector efficiency is given by the following equation: 


pe ee ee (16) 
Qo Gs 


Also, and according to EN12975 [15] for FPCs, the efficiency curve can be written as: 


(AT) CAT)? 
7 Mo ~~ A 2 G (17) 
where: 
1): Optic efficiency 
a, and a, : are the coefficient of the losses according to EN12975 standard. 
G : Global solar radiation (Wm”). 
AT : The temperature difference (K). 


3. METHODS AND MATERIALS 
3.1. Global Algorithm 

Because of the complexity structure and multiplicity of equations of our model, the method of 
Runge-Kutta 4" order [20] was adopted for solving numerically the following systeme of equation: 
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(18) 


The Figure 1 shows the global algorithm following to run and solve numericaly the systeme (18). 
a. Block I: Calculation of the solar radiation 
In this step, we adapted the model of Perrin de Brichambaut [16, 17] as mentioned previously. In the 
beginning we input the parameters of Longitude, latitude, day number and number of the month for 
calculating declination angle and hour angle, and consequently the incidence angle (1) and the high of sun (h) 
are calculated, after this step we input the coefficients exposed in table.2, in the end of this part we have been 
calculated the global solar radiation. 
b. Block II: Calculation of The thermal performance 
Inputting the physics properties (Specific heat at constant pressure, wind speed, thermal 
conductivity and volume density) of water, the geometric parameters exposed of the system. After this step of 


this program starts to calculate the all type heat coefficients (h,,h.,,h.4 ) by adapting the formulas exposed 


cv? 
previously. 
In the end of this step we have adapted the numerical method in order to resolve the system of differential 


equations and consequently the thermal performances (1; Q, ) are calculated. 
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Inputs 
Tam; Wind speed 


Physics properties 
Geometric parameters 


Mass Flow CSS ET A a ee T 1 
Inlet fluid temperature Bloc I: Calculation of the solar radiation by using model of 
P. Brichanbaut [16]-[17]-[18] 





Inputs: 
Day number ; longitude ; 


| Latitude etc 
i Solar Model i 
i i Calculation of the global radiation intensity 
Calculations of all heat coefficients 


Numerical Solution l i Bloc II: Run and solve the system of 
RK-4th i — | eauation include all calculation | 

! | | 

| 








Figure 1. Global algorithm of solving the systemof equation 


3.2. Experimental Set-up 

The experimental apparatus is located on the terrace roof of the Laboratory of Energy of the Higher 
School of Technical Education, Rabat (Latitude: 34 ° O' 47" N; Longitude: 6 °49'57" 0; Altitude 46 m). 
Indirect SWH system with forced circulation consists of 2 m” FPC, a 300L hot water tank heated only by 
solar energy and with an automated sub-system that controlled hot water. In this experiment, the tank heats 
continuously in three days in May and December, without emptied the tank. We preferred to collect data 
every fifteen minutes during the day from Sunrise to sunset. 

Solar water heating system consists of a hot water storage tank (Figure.2.b), control unit, pump 
station (Figure.2.c), and flat plate collector (Figure.2.a). The FPC employed in this study was south facing 
and inclined at 40°. The hot water tank was installed inside the laboratory while the FPC was set up in the 
building plant room. The solar circuits consisted of two types of pipes. The first one is a 4.7 m of copper pipe 
type which is of 18 mm diameter (inside) and 22 mm diameter (outside) that comes from the tank followed 
by a 26 m of caoutchouc pipe type to FPC with 19 mm diameter (inside) and 23 mm diameter (outside). In 
order to minimize the losses of heat, we opted for using Armaflex with 26 mm of diameter (outside) set up as 
insulation for all pipe fittings. The solar circuit pipe length supply and return were 4 m of copper and 25 m of 
caoutchouc. The collecting system consisted of Dietrisol ECO 2.1 flat plate collector with a gross area of 
2.06 m? and aperture area of 1, 91 m°. The solar collector constitutes a box which ensures an effective long- 
term protection, consisting of a framework, an aluminum bottom and successively equipped with the single 
transparent cover glass of 3.2 mm of thickness and an air gap which was located between the cover and the 
absorber which has a 1.2 cm of thickness. The absorber has a geometry constitutes copper pipes welded with 
an aluminum foil where the fluid can be circulated by means of the solar circulator, and an insulator (mineral 
Wool) a 40 mm thickness, the sizes of the FPC is illustrated in table 1. 


Table 1. Specifications of the Different Components of the Solar Collector Used in the Present Study 


Gross dimensions : Width (m) 1,147 Depth (m) 0,087 
Length (m) 1,753 
Transparent cover : Length (m) 1,134 Mass (kg) 16,4 
Material Height (m) 1,740 
Soda lime silicate glass Thickness (m) 0,0032 
Air gap (m) 0,012 Mass air gap 0,027 
(kg) 
Absorber : Length (m) 1,100 
Material Height (m) 1,706 Mass (kg) 10,026 
Aluminum plate Coating: copper Thickness (m) 0,004 
Insulation : Thickness (m) 0,040 
Material With (m) 1,717 Mass (kg) 3,056 
Mineral wool Height (m) 1,1145 
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The working fluids used in solar collectors are mixtures of 30% glycol/water. The collector had 
maximum operating and stagnation temperatures of 164°C and 213°C respectively, a maximum operating 
pressure of 10 bars. The stainless steel hot water tank (model BSC 300) was 1767 mm high with a diameter 
of 601 mm and an operating pressure of 10 bars. The tank had a heating coil with a surface area of 1.2 m2 a 
fluid content of 8, 9 L and a rating of 21. An automated hot water dispensing unit was designed and 
incorporated into the SWHS. The unit includes a programmable logic controller (PLC), contactors, relays, 
electrical fittings, and impulse flow meters. PLC was used to control circulation of the Solar fluid between 
the tank and the PFC, that is when the difference in temperature between the hot water in the bottom of the 
tank and the PFC is superior to 5°C the PLC allows the circulation, while when it is inferior to 3 °C the PLC 
stops the circulation. A pulse flow meter (1 pulse per liter) was used to count the number of liters of Solar 
Fluid from the PFC. 





Figure 2. The test bench, (a) solar collector, (b) storage tank and (c) pump type ST 15/6 ECO, used in the 
experimental measurement 


The SWH is equipped with a RESOL DeltaSol (Figure.3.a) M solar controller which has relay 
inputs to control the operation of the solar pump station. It also has temperature sensor inputs onto which 
PTIOOO platinum resistance temperature sensors are connected to measure water and solar fluid 
temperatures (Figure 1). The volumetric flow rate of the solar fluid is measured using RESOL V 40-06 
impulse flow meters which react at 1 | per pulse. RESOL DL2 data loggers are used to store data every 
fifteen minutes from the RESOL DeltaSol M solar controllers via RESOL VBus cables. DL2 data loggers are 
equipped with a secure digital (SD) drive and a local area network (LAN) port for direct connection to a 
personal computer (PC). Data from the loggers are extracted using a Web browser or an SD card. Global 
solar radiation on the collector's surface is measured using CSIO Solar cell type E (Figure 3.b). 


me ere PS) Tee 





(b) 


Figure 3. the RESOL data acquisition system (a) and cell Type RESOL CS10 (b) 


Experimental and Modeling Dynamic Study of the Indirect Solar Water Heater... (Ouhammou Badr) 


92 o ISSN: 2252-8814 


4. RESULTS AND ANALYSIS 

The main objective of the experimental test is to valid the theory model. In this paragraph, many 
capabilities of our model will be exposed for a sunny and cloudy days. The characteristics of the site and the 
parameters used in the simulation are given in Table 2. The curves exposed in Figure.4 shows the shape of 
the instantaneous efficiency obtained by the theoretical model, the simple experiments and also by the 
EN12975 standard. 
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Figure 4. Collector efficiency versus AT/G. 


The difference between the experimental and the theoretical results is 1,2 % and the same thing 
between the theoretical curve and this one obtained by EN12975 standard. It’s noted that the fluid outgoing 
of the tank and entered the collector has a temperature equal to the first temperature in the beginning; 1.e. we 
always keep the same difference temperatures between the output and input. The difference between both 
curves (experimental and the theoretical) is observed in the high values of AT/G, this means that the present 
system (SWH) may work in high temperature. 

By using the quadratic function, the parameters of efficiencies curves illustrate in equation 23 are exposed in 
Table 2 in comparison with those given by constructor. 


Table 2. Experiments and Simulation results of parameters efficiencies according to EN12975 standard 


Parameters Experiment Simulation the constructor of 
collector 

No 0,774 0,779 0,794 

ay 4.1 3,94 4,31 

ap 0,089 0,09 0,0120 


On the other side, the variation of the outlet fluid temperature during a typical sunny day is 
presented in Figure.5; as it’s shown in the Figureure, the outlet temperature increases at the beginning, passed 
by a maximum and decreases towards the end of the day. 
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Figure 5. Variation of the outlet temperature of the working fluid for a typical sunny day 
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This is due to the evolution of solar intensity flow during the day. The temperature of the fluid 
achieves in maximum about 44.4°C by simulation and 47 °C experimentally. On the other hand, the 
fluctuations in the temperature of the absorber, as it is illustrated in Figure.6, they are caused by the changes 
of the ambient temperatures, glass cover transparent which is caused by the cloudy passages. The percentage 
error concerning the temperature of the absorber is 4 % between the simulation and experimental results. 


TOC) 





500 600 700 #00 $00 1000 1100 


Time (minh 
Figure 6. The variation of the temperature of the absorber of a typical sunny day. 
On the other side, the Figure 7 shows the difference between changes of the useful energy extracted 


from the system. Not only simulation results tend to those of the experiments, but the value of the useful 
energy gain was almost closely predicted by the theoretical model 





(a) (b) 


Figure 7. (a) Solar radiations and (b) Useful energy for a typical sunny day 


As it’s shown in Figure.7 the numerical model value is always under predicts value. Perhaps, among 
the reasons, is that the theoretical model takes the ambient temperature and wind speed as constant 
parameters in the calculation (constant inputs) on the contrary the experiment that has the fluctuation in these 
parameters. The solar intensity obtained experimentally is higher than the one obtained by simulation; this is 
due to the higher value of the Albedo that caused by the existence of the walls and a lot of equipment 
surround our system. 

The performance of FPCs for a cloudy day is shown in Figureures 8-10. The evident point is the 
fluctuation of the instantaneous useful energy gain (Figureure.8.b), that caused by the strong changes in solar 
radiations as shown in the Figureure.8.a. 
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Figure 8. Solar radiations (a) and Useful energy (b) for a typical cloudy day 


The retarded response is very clear in the afternoon. The useful energy and incident radiation 
(experimentally) has a maximum estimated about 1185 W and 1075 W/m’ respectively. The numerical model 
is not well estimated of the useful energy gain extracted from the system for the typical cloudy day. Also, the 
outlet temperature has affected in the same way by the fluctuations of the solar intensity as shown in 
Figure.9. 
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Figure 9. Variation of the outlet temperature of the working fluid for a typical cloudy day. 


Nevertheless, the outlet fluid temperature of the collector was estimated with small discrepancies by 
simulation, 1.e., the maximum experimental value is 42,8 °C when the simulation is 40 °C. It’s necessary to 
note that, when the solar radiations are quickly and unexpectedly increase or decrease, the fluid temperature 
doesn’t change immediately, and according to the changes unexpected of solar radiations, the experimental 
efficiency curve is very different from the simulation as illustrated in Figure. 10. 
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Figure 10. The transient efficiency of the collector for a typical cloudy day 


When the solar radiation is not important, the measurements efficiency is higher than that obtained 
by simulation, the reason responsible for this increases value is because the difference between the inlet and 
outlet experimental temperatures, but for the simulation, as previously mentioned, the parameters inputs are 
constants. In the middle of the day, when the solar intensity increase as shown in Figure.10, the efficiency 
curve drops completely in a different way from that solar radiation. 

The reasons for those deviations between our model and experimental results as illustrate in previous 
Figureures are the following: 
Weather conditions, which were unstable during the experiences. 
The simulation uses all parameters as a constant input. 
The fluctuation in the measurements. 
The distance between collector and storage tank is far. 
The effect of the thermal inertia, i.e. it’s impossible to achieve initial temperatures very close to the 
modeled ones. 
The thermal and optic losses during the experience represent 11 % and 20-- % respectively, for example 
for the cloudy day , the incident energy is 16920 KJ , the losses are 1862 KJ and 3384 KJ respectively , 
for the sunny day the energy incident is 23400 KJ ,the losses are 4680 KJ and 2574 KJ. 

The measuring equipment may not be as accurate as it should be or as indicated by their 
manufacturers, 1.e. the measurements of the solar radiation are carried out with a relative error of 10 % and 3 
% for the other different temperatures sensor. 
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5. CONCLUSION 

The objective of this study was to model an indirect solar water heater with forced circulation and 
validate the model through simple experimental. The simulation is used on this model was shown to agree 
well with the experimental results such as: 

a. The percentage error of outlet water temperature is 2 % for the sunny days and 5 % for the cloudy days. 

b. The percentage error of useful energy is 4 % and 9 % for the sunny and cloudy days respectively. 

c. The thermal efficiency curve of FPCs follows the same shape of EN12975 standard, which the 
coefficients a, and a are equal to those of FPCs (type SOL 200) and the percentage error between them is 
12%: 

This means that the model developed in the current work can be represented the real system and can 
be used to estimate its dynamic behavior. 
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